Search for: All records

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

Structure of the inner core is often measured through traveltime differences between waves that enter the inner core (PKPdf) and waves that travel through the outer core only (PKPab and PKPbc). Here we extend the method to converted waves PKSdf and SKPdf and compare results to PKP wave measurements. PKSdf and SKPdf have a very similar path to PKPdf and if velocity variations are present in the inner core, all three wave types should experience them equally. Since traveltime deviations can be due to velocity changes (either isotropic or anisotropy) as well as wave path deviations born from heterogeneity, we simultaneously investigate wave path directions and traveltimes of PKP, SKP and PKS waves for several source-array combinations. One of the path geometries is the anomalous polar corridor from South Sandwich to Alaska, which has strong traveltimes anomalies for PKPdf relative to more normal equatorial path geometries. Here we use array methods and determine slowness, traveltime and backazimuth deviations and compare them to synthetic data. We find that path deviations from theoretical values are present in all wave types and paths, but also that large scatter exists. Although some of the path deviations can be explained by mislocation vectors and crustal variations, our measurements argue that mantle structure has to be considered when investigating inner core anisotropy. Our polar path data show similar traveltime residuals as previously published, but we also find slowness residuals for this path. Interestingly, SKPdf and PKSdf for the South Sandwich to Alaska path show traveltime residuals that are partly opposite to those for PKPdf, possibly due to an interaction with a localized ultra-low velocity zone where waves enter the core.

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (Sdiff) are potentially strongly influenced by lowermost mantle anisotropy. Sdiff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. Sdiff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SHdiff–SVdiff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of Sdiff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which Sdiff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on Sdiff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SHdiff waves. In particular, we investigate how SHdiff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to Sdiff splitting that can be clearly distinguished from other waveform effects. These results enable us to lay out a strategy for the analysis of Sdiff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of Sdiff and the influence of source-side anisotropy. We demonstrate our Sdiff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained Sdiff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

We compile and make publicly available a global digital database of body wave observations of seismic anisotropy in the D′′ layer, grouped using the method used to analyze deep mantle anisotropy. Using this database, we examine the global distribution of seismic anisotropy in the D′′ layer, evaluating the question of whether seismic anisotropy is more likely to be located at the edges of the two large‐low velocity provinces (LLVPs) in Earth's mantle than elsewhere. We show that this hypothesis lacks statistical justification if we consider previously observed lowermost mantle anisotropy, although there are multiple factors that are difficult to account for quantitatively. One such factor is the global lowermost mantle ray coverage for different phases that are commonly used to detect deep mantle anisotropy in shear wave splitting studies. We find that the global ray coverage of the relevant seismic phases is highly uneven, with LLVP edges and their interiors less well‐sampled than the global average.

Much of our knowledge on deep Earth structure is based on detailed analyses of seismic waveforms that often have small amplitude arrivals on seismograms; therefore, stacking is essential to obtain reliable signals above the noise level. We present a new iterative stacking scheme that incorporates Historical Interstation Pattern Referencing (HIPR) to improve data quality assessment. HIPR involves comparing travel‐time and data quality measurements between every station for every recorded event to establish historical patterns, which are then compared to individual measurements. Weights are determined based on the individual interstation measurement differences and their similarity to historical averages, and these weights are then used in our stacking algorithm. This approach not only refines the stacks made from high‐quality data but also allows some lower‐quality events that may have been dismissed with more traditional stacking approaches to contribute to our study. Our HIPR‐based stacking routine is illustrated through an application to core‐reflected PcP phases recorded by the Transantarctic Mountains Northern Network to investigate ultra‐low velocity zones (ULVZs). We focus on ULVZ structure to the east of New Zealand because this region is well‐sampled by our data set and also coincides with the boundary of the Pacific Large Low Shear Velocity Province (LLSVP), thereby allowing us to further assess possible ULVZ‐LLSVP relationships. The HIPR‐refined stacks display strong ULVZ evidence, and associated synthetic modeling suggests that the ULVZs in this region are likely associated with compositionally distinct material that has perhaps been swept by mantle convection currents to accumulate along the LLSVP boundary.

Shear‐wave splitting measurements are commonly used to resolve seismic anisotropy in both the upper and lowermost mantle. Typically, such techniques are applied to SmKS phases that have reflected (m‐1) times off the underside of the core‐mantle boundary before being recorded. Practical constraints for shear‐wave splitting studies include the limited number of suitable phases as well as the large fraction of available data discarded because of poor signal‐to‐noise ratios (SNRs) or large measurement uncertainties. Array techniques such as beamforming are commonly used in observational seismology to enhance SNRs, but have not been applied before to improve SmKS signal strength and coherency for shear wave splitting studies. Here, we investigate how a beamforming methodology, based on slowness and backazimuth vespagrams to determine the most coherent incoming wave direction, can improve shear‐wave splitting measurement confidence intervals. Through the analysis of real and synthetic seismograms, we show that (a) the splitting measurements obtained from the beamformed seismograms (beams) reflect an average of the single‐station splitting parameters that contribute to the beam; (b) the beams have (on average) more than twice as large SNRs than the single‐station seismograms that contribute to the beam; (c) the increased SNRs allow the reliable measurement of shear wave splitting parameters from beams down to average single‐station SNRs of 1.3. Beamforming may thus be helpful to more reliably measure splitting due to upper mantle anisotropy. Moreover, we show that beamforming holds potential to greatly improve detection of lowermost mantle anisotropy by demonstrating differential SKS–SKKS splitting analysis using beamformed USArray data.

The carbon and water cycles in the Earth's interior are linked to key planetary processes, such as mantle melting, degassing, chemical differentiation, and advection. However, the role of water in the carbon exchange between the mantle and core is not well known. Here, we show experimental results of a reaction between Fe₃C and H₂O at pressures and temperatures of the deep mantle and core‐mantle boundary (CMB). The reaction produces diamond, FeO, and FeH_x, suggesting that water can liberate carbon from the core in the form of diamond (“core carbon extraction”) while the core gains hydrogen, if subducted water reaches to the CMB. Therefore, Earth's deep water and carbon cycles can be linked. The extracted core carbon can explain a significant amount of the present‐day mantle carbon. Also, if diamond can be collected by mantle flow in the region, it can result in unusually high seismic‐velocity structures.